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ABSTRACT 


Simulation generated probabilities Bom an Anti-Submarine 
Barrier composed of submarines are examined from a statisti- 
cal viewpoint. A probabilistic model which is not generally 
considered to be applicable to this case is demonstrated 
to be Statistically Supported. The applicability of the 
model is justified probabilistically. A statistical esti- 
mating relationship is then developed to estimate the sole 


input parameter from submarine FIGURE-OF~MERIT. 
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I. INTRODUCTION 


A. THE TACTICAL PROBLEM 

The encounter treated is that of an Anti-Submarine 
Warfare barrier consisting of submarines. Each barrier 
submarine is confined to an operating area, hereafter called 
a cell, due to the requirement to remain in a safe haven. 
The cell is part of a larger barrier across a restricted 
passage of widen D. The transiting submarine is obliged 
to pass through the passage, and hence the barrier. The 
barrier has uniform cells and may be composed of several 
cells in depth. All detections are made with passive SONAR. 
No external intelligence of the position of either opposing 
force is allowed. All detections are made by the submarines 
and all kills are made without the aid of external forces. 
Transiting submarines are alone, and the time between tran- 
sits is sufficiently long to prevent there being any effect 
of one transit upon another. If cells are arranged in se- 
quence along the transit track, they are sufficiently sepa- 
rated to prevent there being any effect of one engagement 
upon another (This assumption is relaxed for purposes of 
investigation in Section II.C). The barrier is sufficiently 
wide (D is sufficiently large) to make any special effects 
from the cells on the ends negligible. 

By virtue of the above assumptions, the encounter may 


be studied as a single cell with a single transiting submarine, 
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and a Single barrier submarine. The barrier in total may 
then be constructed probabilistically from the building 
block of one cell. Figure (1) is a geographical diagram 
of one cell in the barrier, with pertinent variables 


labelled. 


B. THE SIMULATION 

For the purposes of simulation, the following assump- 
tions were deemed appropriate: 

(i) The transiting submarine chooses a penetration 
point uniformly along the cell front, since the limits of 
the cell are unknown to it. 

(ii) The transiting submarine proceeds at transit speed 
(Vin) until detecting the barrier submarine, at which time 
it slows to quietest speed to evade. 

(iii) The transiting submarine penetrates the cell 
front perpendicularly, and continues on a straight track 
throughout. | 

(iv) No counter measures are employed. 

(v) The barrier submarine proceeds at patrol speed (Vn) 
across its cell until a detection, at which time it commences 
interception tactics. For a short initial period it steers 
the target's bearing, after which time sufficient informa- 
tion is available for a firm track. When firm track is 
established, the barrier submarine proceeds to the intercept 


POint. 
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Figure l 


The Tactical Scenario 
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The simulation was by digital computer. It serves the 
purpose of providing data for analysis only, and is des- 
cribed here to give insight into the likely characteristics 
of the data. The validity of the simulation is not a subject 
of direct analysis. The ability of the data to support | 
reasonable models, and internal consistency were the facets 
of validity tested. The simulation was constructed and 
run by the Center for Naval Analyses and should be con-~ 
sidered separate from the analysis contained herein. 

The simulation is a time-stepped, event-store, Monte- 
Carlo program. Each stochastic Be ecees 1s modeled stochas- 
tically, and each probabilistic event is determined by ran- 
dom number. Input parameters include acoustic character- 
istics of both submarines, sea conditions, initial veloci- 
ties of both submarines and weapons load of both submarines. 
sea conditions were not varied for any simulation runs. 

Two hundred and fifty replications of each submarine 
pairing were conducted. Taking any particular binary value 
output by the simulation to be a bernoulli random variable 
(if.e.: target killed or not killed), the variance of any 


probability estimator is 


Var eD _ on (al = io), 
ri 


This quantity is maximized for a p of .5, so the maximum 


possible variance of such output is .OOl. 


1,6) 





With two hundred and fifty replications the Central 
Limit Theorem will be assumed to hold, and hence the 
estimators will be assumed to be distributed as gaussian 


random variables. 


C. THE ANALYTICAL PROBLEM 

The original analysis of the tactical problem was 
cOnducted by simulation. Seven distinct barrier submarine 
hull and sonar combinations and four distinct transiter 
hull and sonar combinations were included. 

Various quantities were compiled and computed from the 
simulation run, most for the purpose of checking consistency. 
The sole parameters finally desired from this simulation 
were the probabilities of kill for the barrier and the 
transiting submarines. 

The number of distinct simulations is expressed as the 
product of (i) the number of cell widths sufficiently ease 
to allow interpolation and sufficiently numerous to avoid 
extrapolation; (ii) the number of opposing velocity pairs; 
(141) the product of the numbers of distinct opposing 
submarines. The essence of the problem is thus to reduce 
the information contained in such output to a tractable 
analytical expression. The purpose of this reduction is 
four-fold. First, the information must be made convenient 
for a decision maker. Second, the information should be 
expressed as a function of tactical parameters and 


characteristics to allow inspection of the dynamics of the 
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problem and to assess validity. Third, the information 
should be expressed as a function of input parameters and 
decision variables to enable the comparison of alternative 
courses of action. Fourth, an acceptable analytical model 
can either reduce or eliminate the need for further 
Simulation computer time, should additional range of input 


variables be desired. 
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II. THE SEARCH MODEL 


The following discussion treats the searcher as moving 
and the target as stationary, in order to be in consonance 
with discussions in the literature. Clearly, the inter- 
changeability of searcher and target may be used to achieve 
a desired result. 

A common model in search theory is the Random Search 
Model, now a classic in the field [see ref. 1]. The 
assumptions for’this search law to hold are listed in 
ref. l as 

(1) The target (here the searcher) is uniformly 
distributed within A. 

(ii) The observer's (here "the target's") path is 
random in A in the sense of as having its different (not 
too near) portions placed independently of one another 
in A. 

(iii) On any portion of the path which is small 
relatively to the total length of path but decidedly larger 
than the range of possible detection, the observer always 
detects the target (here "thetarget is always detected by 
the peechen™) within the lateral range W/2 on either side 
of the path and never beyond. This last assumption can be 
relaxed to the extent that a number W (search sweep width) 
may be defined by integrating over the lateral range curve 


P(x) where 
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P(x) = P(detect a target passing 
at lateral range x) 


A searcher eakeh lateral range curve P(x) will detect as 
many targets in the entire plane as one with a definite 
range law detector such as indicated by assumption Cini). 
as long as the targets are distributed uniformly, and as 
long as 


co 


f P(x) dx = W 


a= OO 


where W is the diameter of the definite Search law detector. 


Under the Random Search Law, 


P(detection) = 1 = eo WR/A 
W = sweep width 
£ = track length 
A = area containing target 


If the searcher transits through A then 


tI rt 
I 
O}H 


C = width of A 


thus 


_ We 


ll 
} 


P (detection) 
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For tne purposes of this analysis, the assumpticns of 
the Random Search Law hold, with the exception of (ii). 

The searcher (target in this context) is not searching 
along a random track, but 1s actually proceeding directly 
along an effectively zig-zag non-overlapping track relative 
to the target. 

A-second model is developed in ref. 1 which may also 
be applied. This model is a barrier model. It will be 
referred to here as the Channel Patrol Model. 

Suppose a searcher patrols at speed v across a channel 
of width c. Then one patrol (here a channel crossing) is 
completed in = hours. The target is approaching at speed 
u straight down the channel. During one patrol of the 
searcher he travels = miles. Relative to the target, 
the searcher is traveling a Zig-zag path at relative 


speed w, where 


(see Figure 2). 
Thus the total area within which the target may be 


found in one patrol is 


A = (height)(base) = (u =) C 
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Figure 2 


The Channel Patrol Model 
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The total relative area swept by the Searcher in one 


patrol is 


(velocity) (time)(sweep width) 


to 
ll 


(w) (2) W 


If the target is located within this relative area 


uniformly, then the probability of detection is 


u (>) © 


Noting that 
waJfurtv® = uf + (2)? 


then 


P(det) = B/1 + (4)? 


Q| = 


stag 7k << 1 then approximately 


-W 
qk 
P(det) = l-e 


This approximation involves some reduction in the 


probability of kill, since 


W 
“kG 


l =e < k 2 
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This may be justified in that the constrictions of a cell 
or channel will deny the barrier submarine the portion of 
his sensor's range which is beyond the limits of his cell. 
This phenomenon will be called clipping, and will be 
discussed in Appendix B. 

If the probability of kill given detection is p, for 


all targets then 


. W 
“Pak 
P(kill) =l-e 
k= 1 random search 
model 
k = f1 + (Zé channel patrol 
model 


The quantities DPD and k are statistically nuisance parameters, 
and may be suppressed by defining a new quantity, kill 
sweep width, and also denoted W, so that 
(Gl) P(kill) =l-e 
W = kill sweep width 
The advantage of this representation is that p is no 


longer needed, provided that kill sweep width can be 


determined, as will be the case. 
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B. MODEL VERIFICATION 

Rather than questioning the reasonableness of the 
assumptions, or the applicability of the model to the 
tactical prebilen® an alternative approach rene mms b tile. 
Let us suppose the applicability of the model and then 
test the Supposition Statistically. 

It is natural to suppose that W (sweep width) is a 
function of at least (a) speed of both submarines and 
(b) quietness of both submarines. It is explicitly not 
a function of C (cell width). The obvious test of the 
Random Search Model is to invert the formula, injecting 
known C and monte-carlo estimates of P(kill) to obtain 
W estimates, and then test the variability of w with 


respect to C. Now, 


therefore 


W= -C Log, (1 - Py) 


Both a Wr and Wp are found, which are those associated 
with the transiter and barrier respectively. 

From considerations of the data, not included here, it 
is immediately obvious that the variation in W with respect 


to C is very small. This impression was verified statis- 


tically using analysis of variance. It was found that the 
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variations in W were justifiably attributed to statistical 
fluctuation in the monte-carlo estimation of P(kill), and 
not to any systematic model inconsistencies. The detection/ 
kill model ‘ip thus supported statistically (see statistical 
appendix) in the sense that it fails to be contradicted. 

A traditional approach to the simulation data used here 
probably involve the estimation of W from technological 
considerations, and a comparison of this computed value and 
the attendant probabilities to the results of the simulation. 
The approach taken herein is to take the simulation results 
and derive the implied values of W. Validation of the model 
would then involve consistency within the node The main 
reason for this approach is that Seals kill probabilities 
are desired, whereas typicai estimation of w is for Sséarch 
sweep width. Thus inversion of the procedure allows for 
elimination of the necessity of analytically solving the 
kinematical problem of kill, counter-detection and evasion, 
which is the primary source of difficulty in the traditional 
approach. Additionally, should the predicted sweep width be 
sufficiently different from an experimentally derived Ae 
then the classical approach offers little recourse except 
to discard the predicted sweep width, the experimental data 


or the model. 
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C. IMPLICATIONS OF THE MODEL 

There are very important tactical considerations which 
might argue for various barrier configurations for 
submarines. Two building block configurations are here 
examined to determine their probabilistic advisability. 
The first is the "stacked" barrier, with submarines deployed 
in depth. The second is the "packed" barrier with 
submarines deployed along a line (see Fig. 3). 

Now, suppose that there is some larger barrier of 
width D in which these cells are placed. Alternative I is 
to place one packed barrier, alternative II is to stack k 
less dense barriers with the same number of submarines, 


across the channel D. 


Let 
n = number of submarines in the barrier 
P, = P(kill a transiter under I) 
Po = P(kill a transiter under ITI) 
gq, =1l- Py 
G5 1 - Po 


Now, cell width for any single line barrier is related 


to n and D by 


SO 


ay 





ies ae 


eC 
Alternative I 
The Stacked Barrier 


Lim 


eta? mess 
Alternative II 
The Packed Barrier 


Figure 3 
The Packed and Stacked Barriers 
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and assuming independence among stacked barriers 


Wn Wn Wn 
k 


1 - P5 = Io =e e cee A 


where there are k stacked barriers, each with ny submarines 


a = lois .k). Now 


W aa Wn 
- 5 5 = eo 
1 - P5 = 45 = @ = e P = a,=1-p, 
and 
ee ° 


Thus there is nothing probabilistic to distinguish between 
the two barriers under independence. 

Assuming, on the other hand, that there exists some 
Gependence among the k stacked barriers, then either type 
of dependence is possible. The type of dependence results 
from the nature of the opposing forces. There are only 
four possible results from a transit across the barrier. 
They are the following states: 


1 Barrier killed transiter survives 
2 Transiter killed barrier survives 
3 Both submarines killed 

ul No submarines killed 


In the case of 2 and 3, the type of dependence between the 
barrier in question and subsequent barriers is unimportant 
Since the transiter submarine has been killed, and 


thus will not encounter any subsequent barrier. 
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In the case of 1, it is likely that the combat between the 
two submarines would alert the next barrier, or possibly 
Several barriers, and allow for maximum readiness and 
optimum Salinediie® In the case of 4, whatever special 
circumstances allowed the transiter to evade the first 
barrier would be likely to carry over into the second 
barrier. These arguments are expressed analytically as 


follows. Let 


Py; = P(enter state j|leave state i) 


al 
J 


where one state transition occurs upon transiting one 


barrier. Then in these cases 


1,3,4 


}! 
N 
|}! 
Cute 

C4. 

tl 


and 


Pyy > Py, | alae 


Remembering that states 1 and 4 correspond to transiter 
not ‘killed, and states 2 and 3 correspond to transiter 
killed it can be seen that if the first encounter results 
in state 4, then the transiter will be most likely to 
avoid being killed by a state 4 result in subsequent 


barriers. If the first encounter results in state 1, then 


ef 





the transiter is more likely to be killed by a state 2 
result in subsequent barriers. Thus, depending upon the 
likely nature of the first encounter, the decision can be 


made whether to stack or pack barriers. 
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III. PARAMETRIC PREDICTION OF SWEEP WIDTH 


A. DESCRIPTION OF INPUTS 

The applicability of the Random search Model being 
established, the next step is to find a method for predicting 
the value of W in Eq. (1). Since the analysis is intended 
to serve as an aid to decision making, it is desirable to 
use variables in the prediction whose relationship to W is 
simple, and which are under the control of the decision- 
maker, or tactical commander. 

In submarine warfare, an important measure of effective- 
ness of a submarine is its acoustic quietness, called 
FIGURE-OF-MERIT (FOM). FOM measures the self noise of the 
submarine, as well as the signal level of the opponent. 

It combines these quantities into one, expressing the 
quietness of the submarine relative to its opponent. 

The FOM used in this paper is as follows: 


FOM = Le = L, - DT 


3/7 LO Lon’ 1° 


L _/10 
+ 10 Ons) 


2 
! 
Es 
© 
i 
O 
0a 
| 
© 
"oO" 
=) 
© 


Le = Radiated Noise 


radiated machinery noise 


Lem 


1 radiated velocity noise 
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L = Self Noise 


e 
BE ambient self noise 
le. = machinery self noise 
ovis velocity self noise 


DT = Detection Threshold 


Velocities for computation were taken as V the 


B? 
barrier patrol speed, and V,,, the transit speed. After 

detection, velocities changed for both submarines, but no 
account of this was taken, since this represented mainly 


kinematics. 


B. THE ESTIMATING RELATIONSHIP 
The dependent variables to be predicted were Wr and Wp 
(sweep width of transiter and barrier, respectively). 


Independent variables were assumed to be FOM, and FOM 


B Ag 
(FIGURE OF MERIT Barrier and Transiter, respectively). A 
linear relationship was found to be the best predictor of 


Kill Sweep Width, with form: 


= 
tt 


a tb FOMn + qd FOM 


: a,b,d,f,g,h > 0 


= 
il 


f - g FOM, + h FOM, 


B 

Both kill sweep widths were tested, that for the 
barrier (W.) and the transiter (Wn) » using stepwise linear 
regression. In each case, it was found that the most 


Significant variable for predicting sweep width was the 
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FOM of the submarine to whom the W belonged. That is, the 
most Significant variable in predicting Wa was FOM - Boch 
FOM's were found to be significant for both sweep widths. 

The coefficients preceding the variables in Eq. (2) 
were all as might reasonably We expected, except for d, 
the coefficient of FOM, It would seem that the sign of d 
Should be negative, so that improvement in the performance 
of the barrier submarine would degrade the performance of 
the transiter. This would ordinarily be the case, except 
for the implications of the tactical assumptions. Since 
the transiter does not seek combat, all opportunity for 
Such combat arises by way of the barrier detecting the 
transiter. Thus, as the barrier improves in performance, 
more opportunity to attack arises for the transiter. 

In order to ensure adequacy of the regression, two 
additional tests were conducted. For each, the difference 
between the predicted values of Sweep width and the actual 
values, called residuals, were found. These residuals were 
tested using Analysis of Variance. 

The first test concerned the hull types. The question 
tested was whether there was any information associated 
with the Hull/Sonar combinations which was not conveyed by 
means Of FIGURE OF MERIT. The results indicated that there 
was no reason to believe that any additional information 
about the submarine was likely to be of any uSe, beyond 


FIGURE OF MERIT. 
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The second test was whether the patrol speed (Vp) or 
the transit speed (Vin) had any effect beyond that conveyed 
through FOM. The results indicated that the patrol speed 
of the barrier submarine was not significant, whereas the 
transit speed (V,,) was significant. This is due to two 
factors. First, the faster the transiter travels, the less 
time he is exposed to detection. Second, the transiter is 
less susceptible to intercept at higher speed. This 
remaining effect due to Vip is a drawback in the model, 
Since the best possible prediction is not being carried 
out. In view of the excellent fit of the linear regression 
(see the Statistical Appendix) the model seems adequate 
for most needs. It is possible that the additional complica- 


aion of the model or the increased need for model inputs 


would make a more refined model undesirable. 


C. A SECONDARY EFFECT OF THE ANALYSIS 

In the course of the statistical analysis of the 
simulation data, an important result was obtained, above 
and beyond the requirements of data compactness and 
manageability. The result was the detection, from 
statistical considerations, of programming errors in the 
Simulation. Tactical simulations are always large, in a 
computer code sense, and difficult to completely "de-bug". 
In this instance, a logic control variable had been 
included, to allow changing of torpedo P(kill) during a 


run. The program was set up initially using dummy P(kill)'s 
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for all torpedoes. When the final data runs were made, 

the logic control variable was set to inject the correct 
P(kill)'s for barrier types 1 through 4, but was erroneously 
omitted for barrier types 5 through 7. The omission was 

not noticed until the statistical analysis began. As a part 
of the statistical analysis, the residuals (Actual W minus 
predicted W) were plotted for all submarine types. It was 
noticed that barrier types 1 through 4 were all performing 
poorer than predicted, and barrier types 5 through 7 were 
performing better than predicted. This discontinuity was 
not due simply to different Hull/Sonar types, because one 
Hull type and one Sonar were found in both groups, although 
Pot together. When the programmer saw this effect he 
correctiy deduced the error and re-ran the entire set of 
barrier types 5 through 7. The resulting residuals were 
again plotted, and no longer displayed the discontinuity. 

It is unlikely that the error would ever have been discovered 


without the statistical analysis. 
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IV. CONCLUSIONS 


The philosophy pursued herein, that a model is valid 
if it is statistically Supported is a profitable philosophy 
Only insofar as it aids analysis. If a model is reasonable 
from an analytical point of view and also supportable 
statistically, then the best solution has been arrived at. 
It is the case here, as in many instances, that the model 
chosen was not uniquely applicable, nor were its requisite 
assumptions completely satisfied. Hopefully, any doubts 
about applicability were dispelled by the compelling 
Statistical evidence. 

The method of inversion of results to obtain unknown 
inputs, used here to determine W, is also somewhat 
unorthodox. The resultant estimates of inputs are naturally 
only as good as the original results. ASSuming an acceptable 
level of validity, these inverted quantities may be compared 
to engineering or similar estimates, probably to the benefit 
of both the engineering and tactical analyst. Perfectly 
valid and useful models are no doubt frequently rejected 
because the improper estimation of input quantities causes 
inexplicable differences between actual and predicted 
results. It seems reaSonable, however, that statistically 
verifiable consistency in quantities inverted by way of a 
model would argue in favor of the model, and cause 
re-examination of the engineering estimate in case of 


disagreement. Conversely, models which are supposed to be 
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invalid due to assumptions being clearly violated but which 
are found to be statistically supportable can often be 
extended to include the case in question by further 
considerations of the problem. The phenomenon of clipping 


developed here is an example of this potential. 
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APPENDIX A 


STATISTICAL APPENDIX 


1. Validation of the Model 


The formula 


i) = Gea) = ol ae 


= 


kill sweep width 
cell width 


= 


was inverted to yield 
W= - C ln (1 = p) 


Which then becomes a maximum likelihood estimator of W. 


The hypotheses were 


W 
Hy? p=i1- a : 
Hy: Ho is false 
or equivalently 
Hy? W=-=- c¢ In (1 - p) 
Hy: Hy is false 
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The eet was by analysis of variance using a Completely 
Randomized Factorial design with treatments of barrier type, 
transiter type and cell width at seven, four and three 
levels, Fespect ely: Submarine types were taken as 
fixed factors, cell width as a random factor. Two separate 
tests were conducted, one for Wp and one for Wr (barrier 
and transiter kill sweep widths, respectively). F ratios 
and significances were as tabulated in Tables I and II. 
Based on the significance of the cell width treatment, the 


null hypothesis was accepted at Significance level of .05. 
2. Parametric Prediction of Sweep Width 


The estimating relationships 


T atb FOMn ar (al FOM, 


= 
Nl 


+ h FOM 


= 
il 


f - g FOM 


Ab B 


were tested by linear regression. The test statistics for 


Wm were: 
F = 8.9 
R- x 
= .787 
n = 168 


The regression was thus Significant at the level of .05. 
The t¢ statistics associated with the constants b and d 


were 18.714 and 2.98 respectively, both significant at .05. 


By (h 





TABLE I. ANALYSIS OF VARIANCE OF CELL WIDTH FOR Wr 


Treatment — Degrees of iy Significance 


Freedom Ratio 
Barrier ai) 80.1 a 
Transiter 3 153.9 | a 
C 2 2.015 b 


TABLE II. ANALYSIS OF VARIANCE OF CELL WIDTH FOR Wp 


Treatment Degrees of a vignificance 
Freedom Ratio 

Barrier 6 76.94 a 
Transiter 3 181.621 a 
C 2 1.6515 b 
Notes: 

a << .001 

Da? “se 
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In each regression, data included each of twenty-eight 
submarine pairings at two separate values of Vp and three 
separate values of V,, making a total of 168 data points. 
In each regression there were no other variables ever 
tested as independent variables, thus the significance 
Meveireported is valid. 

Upon completion of each regression, the residuals were 
tested for any variability remaining which could be removed. 
This was done by considering two analyses of variance of 
each set of residuals. 

The first post-regression analysis of variance used 


barrier and transiter hull/sonar type as treatments at 


7 and 4 levels respectively. The hypotheses were 


H,: Hull/Sonar type has no effect on the residual 


Ho false 


Test design was Randomized Block Factorial Design with 
a block consisting of the six distinct velocity pairings. 
The two treatments were considered as fixed. F ratios and 
significances were as tabulated in Tables III and IV. Based 
on the resulting significances, the null hypothesis was 
accepted. The overall significance was .05. This implies 


an individual significance of .013, where 


P(No test value outside critical region) = are 
thus, (1-4)) = .95 
0 = 1- (.95)*= .013 


oe, 





TABLES IT. 


Treatment 


Barrier 


Transiter 


TABLE IV. 


Treatment 


Barrier 


Transiter 


ANALYSIS OF VARIANCE OF RESIDUALS QF Wm 
(SUBMARINE-TYPE TREATMENTS) 


Degrees of 
Freedom 


6 
3 


ANALYSIS OF VARIANCE OF RESIDUALS OF W 


(SUBMARINE-TYPE TREATMENTS ) 


Degrees of 
Freedom 


a 


) 


3 


40 


Significance 


.8173 
~ 3991 


B 


Significance 


EST 
.0216 





The second post-regression analysis of variance used 
Vp and Vin as treatments at two and three levels respectively. 


The hypotheses were: 


Hy: Vp has no effect on the residuals 


H H,~ false 


a 0 
and 


Pies: Vin has no effect on the residuals 
i, : Ho false 

For this experiment, there were four separate tests, 
one for each unit's velocity, and thus two for each 
regression. The likelihood is that one or more of these 
null hypotheses might be true with no implications for the 
others. Experiment design was again Randomized Block 
Factorial with each effect now taken as randon. Fach block 
consisted of 28 distinct submarine pairings. F ratios and 
Significance levels are tabulated in Tables V and VI. 
Based on the resulting significances, the velocity of the 
transiting submarine is significant for the probability of 
its detection and kill by the barrier submarine. This 
result 1s probably due to the resulting reduced time 
available to the barrier submarine to detect, and to the 
increased difficulty in closing with and killing given 


detection. 
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TABLE V. ANALYSIS OF VARIANCE OF RESIDUALS OF W 


(VELOCITY TREATMENTS ) 2 
Treatment — Degrees of i Significance 
Freedom Ratio 
Vp ak Si. 20 
Vin 2 7.4 12 


TABLE VI. ANALYSIS OF VARIANCE OF RESIDUALS OF W 


(VELOCITY TREATMENTS) 2 
Treatment Degrees of F Significance 
Freedom Ratio 
Vp 1 6.1 aL 3 
Vin 2 24.7 ~039 
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APPENDIX B 


THE CLIPPING PHENOMENON 


Suppose that a searcher has Lateral Range Curve p(x) 
(vanishing after R) the probability of detecting a target 


passing at lateral range x. Define 


R 
W=2f p(x) dx 
0 


Then if the target is uniformly distributed across width C, 


centered at the observer, 


ol E> 


P(detect) = R=. 'C 


Now, however, suppose that the curve p(x) vanishes 
outside a range R greater than C, and that the searcher is 
constrained to remain within width C and thus may not 
pursue nor take any action with respect tova target outside 
of Cc. Let the searcher's position from the edge of the 
channel be denoted u (Figure 4). 

In this case, the probability of detecting a Baeest 
with lateral range (closest point of approach) distributed 


uniformly across C, for a searcher at u is 


P(det) = f p(x) = ax + ff p(x) = dx 
, 0 
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No clipping, searcher centered in channel 


p(x) 


u 


Clipping, searcher at u 


Figure 4 


Crip pane 


yy 





denoted hereafter as 
P(det) = 2 {P(u) - P(c-u)} = Sw (u) < dw 


where 


u 
Boy = f- p(x) dx 
0 


Now, if the searcher is in turn uniformly distributed 


in position u across C, the probability of detection is 


le _ i. 
P(det) = § J P(x) - Ple-x) ax = 3W, 
, 


Now, if the searcher is located in the center of the 


channel, 


and so 


- - i C 
P(det) = a W, (5) 
Clearly, if p(x) is decreasing in x for x positive, and 
symmetric then 
C 
W (5) 


C— C 





W 
<¢ 


Thus against a uniform target, probability of detection 


is greatest for a central position, and is diminished by 
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patrolling. Of course, if patrol is not conducted, the 
possibility arises that barrier station location may become 
compromised, so a patrol is clearly necessary. 


Let p(x) be of the form 
p(x) = e *% a> 0 


which is reasonable for passive sonar detection (see Ref. 1), 
but intended only as an illustration. If the searcher is 


at u = = (centered in the channel) 


a, 
P(det) = ew, ($) =e -e ") 


If the searcher is uniformly distributed across the channel 


2C 


le = Mae cK W 
P(det) = ¢ Wa = & (1 sel 2 ye) 
These quantities are graphed in Fig. 5 as a function of 


es along with the quantities 


_W W 
Ps CG a<t 

ail W 

= J] aq? tl 
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The agreement is seen to be excellent among the first three in 
areas of mild clipping, and when there is excessive clipping, 
the Random Search Law (Equation (1)) is somewhat optimistic. 
This contrasts sharply with the obvious inaccuracy of the 
last quantity, the Definite Range Law. It should now be 
fairly clear: that the amount of clipping depends a great 
deal on the shape of p(x), the lateral range curve. 
Lacking any better information, the slight degradation 
implied in the Random Search Law is quite probably suffi- 
cient to a reasonable degree of accuracy. If, as in this 
analysis, the data support the model, then there should be 
no doubt as to its reasonableness. 

Most cases of barrier patrol will probably be found to 
lie between the cases of the searcher centered in’ the 
channel and the searcher distributed uniformly across a 
cell. The problem treated here was within the domain of 
excellent agreement (moderate clipping), and hence the 
law should reasonably hold, since the Lateral Range Curve 
Should have been roughly exponential in x. Thus the 


applicability of the formula 


~ 


Ws 


P(kill) =l-e 


is established. 
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